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■ ABSTRACT 

I Accurate calibration of data is essential for the current generation of CMB experi- 

lO ■ ments. Using data from the Very Small Array (VSA), we describe procedures which 

OO ' will lead to an accuracy of 1 percent or better for experiments such as the VSA and 

\ CBI. Particular attention is paid to the stability of the receiver systems, the quality 

t:;^ . of the site and frequent observations of reference sources. At 30 GHz the careful cor- 

' rection for atmospheric emission and absorption is shown to be essential for achieving 

OO , 1 percent precision. 

■ The sources for which a 1 percent relative flux density calibration was achieved 
^-J I included Cas A, Cyg A, Tau A and NGC7027 and the planets Venus, Jupiter and 

• i-H , Saturn. A flux density, or brightness temperature in the case of the planets, was 

' derived at 33 GHz relative to Jupiter which was adopted as the fundamental calibrator. 

H \ A spectral index at ~ 30 GHz is given for each. 

^ Cas A, Tau A, NGC7027 and Venus were examined for variability. Cas A was 

found to be decreasing at 0.394 ± 0.019 percent per year over the period March 2001 
to August 2004. In the same period Tau A was decreasing at 0.22 ± 0.07 percent 
per year. A survey of the published data showed that the planetary nebula NGC7027 
decreased at 0.16±0.04 percent per year over the period 1967 to 2003. Venus showed an 
insignificant (1.5 ± 1.3 percent) variation with Venusian illumination. The integrated 
polarization of Tau A at 33 GHz was found to be 7.8 ± 0.6 percent at pa = 148° ± 3°. 

Key words: cosmology: observations - cosmic microwave background - techniques: 
interferometric - methods: observational - radio continuum: Solar system - radio 
continuum: ISM 



1 INTRODUCTION 
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temperature scale of each experiment. The cosmological sig- 
nificance of the data is directly dependent on this scale. For 
example, the amplitude of the first peak in the CMB power 
spectrum is proportional to 57b while the ratio of the first to 
the third peak, derived from a mix of experiments probing 
different angular modes I, is also sensitive to Q.b. Further- 
more, the fractional uncertainty in the power spectrum is 
twice that of the fractional uncertainty in the temperature 
scale. Current CMB experiments would greatly benefit from 
a coherent calibration across experiments such as Planck, 
WMAP, CBI, VSA etc. at an accuracy of 1% or better. The 
present investigation provides a set of source intensities cal- 
ibrated to this precision at a frequency of 33 GHz. 

We have chosen the brightest planets and strongest ra- 
dio sources to calibrate the VSA on a frequent (typically 
few hours) basis during CMB observations. The basic cal- 
ibration of all the targets used in the VSA is in terms of 
an assumed brightness temperature of Jupiter (Watson et 
al. 2003; Dickinson et al. 2004), which in turn, is tied to the 
CMB dipole (Hinshaw et al. 2008). 

Inevitably the strongest sources used for the frequent 
calibration of VSA data are sufficiently extended (a few ar- 
cmin) that their total flux densities require correction when 
applied to the ~10 arcmin resolution of the VSA at an ac- 
curacy of better than 1%. Cas A, Tau A and Cyg A fall into 
this category although for Cyg A the effect is ^0.3%. The 
planets do not require such a correction. 

The polarization of the calibrators also needs to be 
taken into account. The most strongly polarized source is 
Tau A. The integrated emission from HII regions and the 
planets are not expected to be polarized. The radial polar- 
ization of the planets would only be of concern if there were 
a significant phase effect in the brightness distribution across 
the planet. 

Time variability is an important consideration when us- 
ing calibrators. Many of the stronger extragalactic sources 
such as quasars have to be rejected for this reason. This is 
an integral part of the present study. 

The spectral index of a calibrator needs to be speci- 
fied if the calibration process is to be useful between exper- 
iments which have different central frequencies and band- 
widths even within the 1 cm wavelength band investigated 
here. In general, spectral index data are taken from the lit- 
erature. 

The main challenge in accurate calibration of ground- 
based observations at centimetric and millimetric frequen- 
cies is provided by the atmosphere. A high dry site is a 
prerequisite. The Izafia observatory of the lAC has a proven 
record of observations in the range 10-33 GHz (Davies et 
al. 1996). An assessment of its properties at these frequen- 
cies will be given in a separate paper (Davies et al., in prep.). 

This paper is arranged as follows. Section 2 describes 
the relevant features of the VSA used in this investigation. 
The philosophy of the approach, including the method of de- 
termining the atmospheric corrections, are given in Section 
3. The results given as flux densities, or brightness tempera- 
tures in the case of the planets, are presented along with the 
adopted spectral indices in Section 4, which also includes a 
discussion of source variability and polarization. The con- 
clusions of this work are summarised in Section 5. 



2 THE VSA AND ITS CALIBRATION SYSTEM 

The advantage of using an interferometer for "point source" 
flux density measurements is that the extended emission 
from the sky, atmosphere and the ground can be largely 
eliminated. The VSA is described in detail by Watson et al. 
(2003), Scott et al. (2003) and Dickinson et al. (2004). A 
summary of its main features is now given. 

2.1 The overall system 

The Very Small Array (VSA) is a 14-element interferometer 
operating in the Ka band (26-36 GHz). It is located on the 
high and dry site at the Teide Observatory, Izafia, Tener- 
ife at an altitude of 2340m. The antenna of each element 
consists of a conical corrugated horn feeding a paraboloidal 
mirror. Each antenna and its associated cryogenically cooled 
receiver system can be placed anywhere on a 4-m x 3-m 
tip-tilt table located in a metal enclosure to minimise any 
emission contribution from the ground. Due to the geome- 
try of the table and the enclosure, the VSA declination is 
restricted to the range —5° to -1-60°. Tracking of a field is 
accomplished by a combination of table tilt and rotation of 
each mirror. 

Since the beginning of CMB observations in September 
2000 the VSA has been operated in both a compact and 
extended configuration. In the compact configuration the 
mirror diameters were 143mm giving a primary beamwidth 
(FWHP) of 4?6. The geometrical arrangement of the anten- 
nas on the table gave interferometer baseline lengths ranging 
from 0.20m to 1.23m allowing angular multipoles of 1 = 150 
to 900 to be probed at the observing frequency of 34.1 GHz. 
The resolution of the array was ~ 17 arcmin depending on 
the field declination and hour angle range observed. 

The extended array, which became operational in Oc- 
tober 2001, used an antenna configuration which covered a 
baseline range of 0.6m to 2.5m thereby extending the range 
of angular multipoles to 1 = 1500 and gave a synthesized 
beamwidth (FWHP) of ~ 11 arcmin. The larger mirror di- 
ameters provided a factor of 1.6 improvement in filling factor 
and an overall factor of ~ 3 in temperature sensitivity (see 
Table [l|. Since September 2005 the VSA has operated in 
a third configuration, the super-extended array, with 60 cm 
aperture mirrors (Genova-Santos et al. 2008). 

The VSA has operated at either 33 or 34.1 GHz with 
an instantaneous bandwidth of 1.5 GHz. With an average 
system temperature of some 35 K the VSA achieves an 
overall instantaneous point source sensitivity of « 36 and 
~ 6 Jy s"^^^ in the compact and extended configurations re- 
spectively; the corresponding temperature sensitivities are 
40 and 15 mK s^'''^. The VSA is linearly polarized in the 
vertical direction and is therefore sensitive to any linear po- 
larization in the calibration sources. 

2.2 Present analysis 

An interferometer array like the VSA requires an accurate 
determination of the position of the antennas to produce 
reliable relative observations at the 1% level aimed at in this 
paper. A maximum likelihood technique has been applied 
to the calibrator source observations in order to constrain 
such telescope parameters as the antenna positions, effective 
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observing frequencies and correlator amplitudes and phase 
shifts (Maisinger et al. 2003). The stability of the VSA is 
such that these calibration observations are required only 
every few hours. 

The overall gain of each antenna is monitored contin- 
uously by means of a noise injection system. A modulated 
noise signal is injected into each antenna via a probe in the 
horn and is measured using phase-sensitive detection after 
the automatic gain control (AGC) stage in the IF system. 
The relative contribution of the constant noise source to the 
total output power from each antenna varies inversely with 
system temperature and thus a correction can be made to 
the overall flux calibration. This system allows account to 
be taken of the gain and variations in total set noise due, for 
example, to atmospheric emission. It potentially provides an 
indication of weather conditions and is a primary indicator 
for flagging weather-affected data in CMB observations. In 
the face of the stringent requirements here of better than 
1% calibration accuracy, this approach had to be modified 
as described in the next section. 



3 THE SOURCE CALIBRATION 
PROGRAMME 

We consider here the methods used to measure relative flux 
densities to an accuracy of 1 percent or better at a frequen- 
cies around 30 GHz. 



3.1 Calibration sources for CMB observations 

The sources available for this study are those bright enough 
to be detected with good signal to noise in say 20 minutes of 
observing time. This limits them to 10 Jy in the two con- 
figurations, which allows the following sources to be used as 
primary calibrators: Tau A, Gas A, Gyg A, Jupiter, Saturn 
and Venus. Observations on these sources of 10-30 minutes 
are interleaved within GMB fields. These observations pro- 
vide the data archive on which the present investigation is 
based. Additional long-track observations covering ±3 hr in 
hour angle were used to quantify elevation-dependent atmo- 
spheric effects. 

3.2 Data reduction procedures 

The basic data reduction of VSA data followed the same 
procedure as described in early VSA work (e.g. Taylor et 
al. 2003; Dickinson et al. 2004). One radio source (or the 
CMB) is calibrated against another (the " calibrator ") by 
converting the units of the 91 raw visibilities from correla- 
tor units to to flux densities via the primary " calibrator " 
for each baseline and then combining them to give a mean 
flux density for the source. The resulting 91 baseline visi- 
bility tracks are checked by eye and interactively searched 
for residual contaminating signals. Due to the large num- 
ber of calibrator observations we were able to be fastidious 
about data quality. Some 20% of the data are removed in this 
way in addition to the more obvious errors due to hardware 
failure and bad weather. The elevation of each source was 
recorded and could be used to make subsequent corrections 
for atmospheric effects. 



3.3 Choice of fundamental reference source 

The flux calibration of VSA data is tied to observations of 
Jupiter which is widely used as a calibrator for GMB obser- 
vations in the range 10-500 GHz. Our early VSA results (see 
e.g. Scott et al. 2003; Taylor et al. 2003) were based on a 
brightness temperature Tj = 152±5 K at 32 GHz, extrap- 
olated to 33/34.1 GHz using 13 = -2.24 ± 0.1, from direct 
measurements (Mason et al. 1999). The VSA calibration is 
now expressed in terms of the WMAP spectrum of Jupiter 
which is tied to the GMB dipole temperature. At 33 GHz, 
the most recent WMAP value is Tj = 146.6±0.75 K (HiU et 
al. 2008)0 The ultimate absolute precision of our adopted 
intensity scale is therefore « 0.5 %. Of course we can, in 
principle, measure relative flux densities/brightness temper- 
atures more accurately than this. 

Accurate determination of relative intensities of sources 
is made from the daily calibration observations and includes 
data from both the compact and the extended VSA configu- 
rations. Longer term flux density or brightness temperature 
variability measurements are possible from the observations 
which extend from September 2000 to October 2004. 

The brightness temperature spectrum of Jupiter is 
shown between 0.1 and 200 GHz in Fig. [1] The broad de- 
pression at around 20 GHz is believed to be due to molecular 
absorption by ammonia. Fig. [1] also shows this part of the 
spectrum in more detail mainly delineated by the accurate 
WMAP points (Page et al. 2003). The best-fit temperature 
spectral index at 33 GHz is /? = -(-0.248±.004. Models for the 
Jupiter spectrum by Gulkis et al. (1974) and Winter (1964) 
give an indication of the expected shape in this frequency 
range. 

3.4 The atmospheric contribution 

The major effect of the atmosphere on the calibration pro- 
cess in the VSA is its variable contribution to the total sys- 
tem noise which is subject to the automatic gain control 
(AGG). The second smaller contribution is the absorption 
of source radiation on its passage through the atmosphere 
to the antenna. Both effects vary with the cosecant of the 
elevation. A detailed discussion of the atmospheric contri- 
bution is given separately by Davies et al. (in prep.). 

The AGC system gives a reduction of the astronomical 
signal relative to the zenith value by a factor, 

/l = [Til 4- Tatm(2)]/[T_R + Tntm(E)] (1) 

= 1/[1 + (T„t™(.)(cosec£ - 1)/{Tr + T„t™(.)))] 

where T_r is the receiver noise plus cosmic background 
and Tatm{z) and Tatm.(E) are the atmospheric emission at the 
zenith and at elevation E respectively. For an average total 
set noise at the zenith of T_r -|- Tatm{z) = 35 K applicable 
to the VSA antennas and Tatm(z) = 7.5 K, we obtain 

fi~l- 0.23(cosec£ - 1) (2) 

^ It is of interest to note that the previous WMAP value (Page 
et al. 2003) is identical to the new one (Hill et al. 2008) but 
with a larger error (146.6 ± 2.0 K) due to beam uncertainties. 
This led the CBI group (Readhead et al. 2004) to use a weighted 
estimate between the Mason et al. (1999) and earlier WMAP 
values, obtaining Tj = 147.3±1.8 K at 32 GHz. 
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Table 1. Specification of the VSA in the compact and extended configurations. 





Compact 


Extended 


DecUnation range 


-5 to +60 


-5 to +60 


Number of antennas (baselines) 


14(91) 


14(91) 


Range of baseline lengths 


0.20 m to 1.23 m 


0.6 m to 2.5 m 


Centre frequencies 


34 GHz 


33, 34 GHz 


Bandwidth of observation 


1.5 GHz 


1.5 GHz 


System temperature, Tsya 


» 35 K 


^ 35 K 


Mirror diameters 


143 mm 


332 mm 


Primary beam (FWHM) 


4.6 


2.0 


Synthesized beam (FWHM) 


» 17 arcmin 


11 arcmin 


Range of angular multipole (1) 


150 to 900 


300 to 1500 


Point source sensitivity 


Si 36 Jy sl/2 


« 6 Jy sl/2 


Temperature sensitivity per beam 


^ 40 mKsi/2 


f» 15 mKsi/2 


Polarization 


linear (vertical) 


linear (vertical) 
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Figure 1. The brightness temperature spectrum of Jupiter. (a) 
The spectrum covering 0.1 to 300 GHz from published data; cy- 
clotron emission dominates the lower frequencies and thermal 
emission the higher frequencies, (b) The spectrum centred around 
30 GHz showing the accurate observations from WMAP (Page et 
al. 2003). Models by Gulkis et al. (1974) and Winter (1964) show 
the effect of ammonia absorption. For the model of Winter (1964), 
we used values of Tn = 168 K and Pn = 2 atm. 



comprising 6.5 K of O2 emission and 1.0 K water vapour 
emission, equivalent to 3mm precipitable water vapour. A 
similar factor relating the intensity of the astronomical sig- 
nal arriving at the antenna from an elevation E to that which 
it would have been at the zenith is 



/2 = [e 
= 1 



-r(B)]/[g--(^)] 

- r(2;)(coseci5 — 1) 



(3) 



For the relevant model atmosphere above Izana the at- 
mospheric temperature is 260 K, so with Tatm{z) ~ 7.5 K, 



f2 = l- 0.028(cosec£ - 1). 



(4) 



The final correction to the source observations for the 
atmosphere, to first order in (coseciS — 1), is 



/ = /i X ,f2 = 1 - 0.26(cosec.B - 1). 



(5) 



Since the VSA is restricted to elevations above 50° , the 
corrections are less than 7% and the approximations in the 
expressions for /, /i and /2 are valid. In each pair of observa- 
tions for a fiux density comparison, each source is corrected 
for elevation as in equation (5) and the source ratio corrected 
to the zenith is obtained. 

A convincing demonstration can be made of this cali- 
bration technique. The quality of the data corrected as de- 
scribed above can be assessed by intercomparing the 3 cali- 
brators Jupiter, Cas A and Cyg A, all observed near transit 
and spread over the period 2002 to 2004. Fig.[2]shows plots of 
the observed flux density of Cyg A and Cas A using Jupiter 
as a calibrator. A similar plot for Tan A is given in Fig. [9] 
below. The daily scatter in the fluxes measured relative to 
Jupiter is seen to be < 1% . The level of thermal noise in a 
typical daily measurement is ^ 0.1 Jy while the daily scat- 
ter is ~ 1 — 2 Jy, thus the variations are dominated by the 
atmosphere. In the case of Cyg A and Cas A, the time dif- 
ference between the source and Jupiter observation can be 
as large as 12hr. For Tan A, where the time difference is 
smaller (< 6hr) over this period, the scatter was < 0.5%. 

Since the main contributor to the scatter is atmospheric 
water vapour emission, the 1% scatter corresponds to 1% 
of set noise within the AGC system, which is 0.35 K. At 
33 GHz this temperature scatter corresponds to 1 mm pwv 
(Danese & Partridge 1989). Some of the implications of this 
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Figure 2. VSA flux density measurements of Cyg A and Cas A 
calibrated by Jupiter over the period 13 March 2001 to 9 August 
2004. The solid lines are least-square fits to the data. 



high level of stability in the atmosphere is included in the 
discussion of individual sources in sections 4.1 to 4.3. 



4 DERIVED SOURCE PARAMETERS 

We now derive the flux density and/or brightness temper- 
ature of our calibration sources at 33.0 GHz relative to 
Jupiter. 

Since our results are most useful to the community 
working in the ~ 10 GHz wide band around 33 GHz, 
we have derived a spectral index for that range from flux 
densities in the literature (plus our new determinations). 
The brightness temperature spectrum of Jupiter has already 
been discussed in Section 3.3. We extrapolated all the VSA 
data to a common frequency of 33 GHz using the spectral 
indices derived from the literature. 

For several of the sources in our study, data are available 
for time spans between 1 and 4 years. These accurate data 
are used to investigate source variability at high precision. 



4.1 Cygnus A 

Cyg A is a radio source associated with a double galaxy at a 
redshift of z = 0.056 containing a central core and hot spots 
at the outer ends of diffuse radio lobes. The overall extent is 
2.1 arcmin and is effectively unresolved in the compact and 
extended VSA arrays. No variability has ever been reported 
for Cyg A at any wavelength. 
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Figure 3. The flux density spectrum of Cyg A from 10 MHz to 
100 GHz; the low frequency turnover is due to Galactic free-free 
absorption. 



Table 2. Cassiopeia A flux density for epochs 1965 and 2000. 



Frequency 


Sl965 


S2000 


Secular decrease 


GHz 


Jy 


Jy 


% per year 


1.405 


2439±50 


1970±50 


0.67 


1.415 


2470±50 


1885±40 


0.67 


1.44 


2328±50 


1793±40 


0.66 


1.44 


2367±20 


1813±55 


0.66 


3.15 


1258±38 


1030±30 


0.62 


3.2 


1279±58 


1008±45 


0.61 


4.08 


1084±26 


860±20 


0.59 


6.66 


684±20 


548±16 


0.57 


8.25 


615±22 


497±18 


0.55 


10.7 


468±0 


388±0 


0.54 


13.49 


394±13 


323±11 


0.53 


14.5 


367±10 


310±9 


0.52 


15.5 


376±18 


309±15 


0.51 


16 


354±11 


292±9 


0.51 


22.28 


285±10 


236±9 


0.49 


32 


224±6 


192±5 


0.47 


33 


211±5 


183±5 


0.47 


86 


115±4 


100±4 


0.41 


87 


109. 4±0 


95.4±0 


0.41 


140 


78.3±7 


69.1±6.2 


0.38 


250 


51.8±5.6 


47.2±5.3 


0.36 



The radio spectrum of the integrated emission of Cyg A 
taken from literature is shown in Fig. [3] There are several 
changes of slope in the spectrum over the frequency range 0.1 
to 100 GHz; one is at ~ 0.5 GHz and another is at ~ 3 GHz. 
The slope of the section 5-100 GHz, which includes 33 GHz, 
is —1.208 ± 0.017. We recommend this as the spectral index 
at 33 GHz. This value may be compared with -1.244±0.014 
given in Table 3 of Baars et al. (1977) at similar frequencies. 

Observations of the calibration sources Cyg A, Cas A, 
and Jupiter were analysed for the period March 14 2001 to 
4 May 2004. The flux density of Cyg A cahbrated by Jupiter 
(assumed to have a brightness temperature of 146.6 K at 33 
GHz) over this period is shown in Fig [51 In this period the 
RA of Jupiter ranged from RA=08h to 12h so observations 
were typically separated by 8-12 hr. The 1% scatter in the 
derived flux density of Cyg A was due to the atmosphere 
(mainly water vapour) over this time range. The mean flux 
density at 33 GHz was found to be 36.4±0.2 Jy. 
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These data can also be used to probe the long-term 
stability of the radio emission from Cyg A and Jupiter at 
33 GHz. A fit to the data shows that the flux density ratio 
of Cyg A and Jupiter increases by 0.043±0.039% per year. 
Accordingly we conclude that the relative radio emission 
changed by ^ 0.1% per year over the period March 2001 to 
May 2004. 



4.2 Cassiopeia A 

Cas A is a 330- year old shell-type SNR of 5 arcmin diameter. 
Baars et al. (1977) give its spectral index between 0.3 and 31 
GHz as a = -0.792±0.007 in 1965 and a = -0.770±0.007 in 
1980. Radio imaging shows multiple hot spots with spectral 
indices between —0.65 and —0.90 (see e.g. Anderson & Rud- 
nick 1996) embedded in extended structure. Steeper spectra 
are associated with features thought to be bowshocks and 
with features outside the main radio ring; flatter spectra are 
found in the ring and in bright features within it which may 
account for its slow but measurable decrease in flux density 
with time at a rate which is thought to vary with frequency. 
These are important considerations when using Cas A as a 
calibrator. 



4-2.1 The Cas A spectrum 

Fig. |4] shows the spectrum of Cas A taken from published 
data over the frequency range 1.0 to 100 GHz with the 
flux densities corrected to the years 1965 and 2000 using 
the recipe given in Section 4.2.3 and listed in Tabled The 
spectrum over this large frequency range shows a flatten- 
ing of the spectral index above 15 GHz not detected in 
earlier data; this flattening is seen in both the 1965 and 
2000 data. In 1965 a = -0.795 ± 0.12 at < 15 GHz and 
a = -0.725 ± 0.013 at u > 15 GHz. For the 2000 data 
Q = -0.767 ± 0.018 at < 15 GHz and a = -0.663 ± 0.012 
at > 15 GHz. An appropriate spectral index at 33 GHz at 
the present epoch would appear to be q = —0.69 ± 0.02. 



4-2.2 VSA measurements of Cas A 

With the resolution of 17 and 11 arcmin for the compact 
and extended forms of the VSA, Cas A will show significant 
resolution effects in a project which aims for an accuracy 
of better than 1%. Accordingly, a correction was applied 
in the analysis procedure which scaled the visibilities on 
each baseline. The 5 GHz map from the VLA was used for 
this purpose. In view of the differing spectral indices within 
Cas A, it was necessary to confirm that the 5 GHz model was 
applicable to our 33 GHz observations. Fig. [S] which plots 
the visibility amplitude of Cas A versus the baseline length 
shows the results of this investigation. It can be seen that 
the normalised visibility after correction with the model is 
unity within 1% across the baseline range 20 to 250 A. The 
compact and extended VSA data sets overlap at a precision 
of better than 0.5% in the baseline range 50-130 A. The cor- 
rection model is thus adequate for our purpose. Fig. [IJb) 
shows the flux density of Cas A measured over the period 
13 March 2001 to 9 August 2004. Each point is the aver- 
age of the flux densities measured relative to Cyg A and 
Jupiter. The secular decrease in the flux density of Cas A is 
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Figure 4. The spectrum of Cas A for the epochs 1965 (crosses) 
and 2000 (filled circles) calculated from published data with the 
secular decreases applied appropriate to each frequency (see Sec- 
tion 4.2.3). 
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Figure 5. The visibility amplitude of Cas A as a function of 
VSA baseline length measured in wavelengths. The crosses show 
the visibilities measured with the extended array; the filled circles 
show the compact array. The open circles are the visibilities cor- 
rected to a 5 GHz observed map of Cas A; this correction brings 
the visibility to within 1% of unity. 



clearly seen. No seasonal effects are evident in the data as 
the 3 sources move from day to night. Moreover the data for 
the compact array (March 2000 to May 2001) are entirely 
consistent with those from the extended array; no step is 
seen in the plot around May 2001. A linear fit to the data 
of Fig. [2][b) gives a flux density for Cas A on the adopted 
Jupiter scale of 182.02±0.07 Jy on January 2001. The linear 
rate of decrease over the period March 2001 to August 2004 
is given by 



dS 
Sdt 



-0.394 ± 0.019% yr" 



(6) 
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Table 3. Secular decrease of Cas A as a function of frequency. 



Frequency 


Epoch 


Decrease 


GHz 




% per year 


0.038 


1955-87 


0.8±0.08 


0.0815 


1949-69 


1.29±0.08 


0.0815 


1949-1989 


0.92±0.16 


0.0815 


1965-1989 


0.63±0.06 


0.1025 


1977-1992 


0.80±0.12 


0.927 


1977-1996 


0.73±0.05 


0.950 


1964-1972 


0.85±0.05 


1.405 


1965-1999 


0.62±0.12 


1.420 


1957-1976 


0.89±0.02 


1.420 


1957-1971 


0.89±0.12 


3.000 


1961-1972 


0.92±0.15 


3.060 


1961-1971 


1.04±0.21 


7.800 


1963-1974 


0.7±0.1 


9.400 


1961-1971 


0.63±0.12 


15.500 


1965-1995 


0.6±0.06 


33.000 


2001-2004 


0.39±0.02 



4-2.3 Secular decrease in flux density of Cas A 

The current estimate of the secular decrease in the flux den- 
sity of Cas A at 33 GHz obtained over a 3.5 year period in the 
present decade is the most accurate of any epoch. Scott et 
al. (1969) at 81.5 MHz covering the two decades 1949-1969 
found a value of 1.29 ± 0.08% yr'^ Using more recent data 
at 81.5 MHz, Hook et al. (1992) found a lower value rate of 
decrease between 1949 and 1989.7 of 0.92 ± 0.16% yr"\ 

We now re-examine the data available relating to the 
variability of Cas A as a function of both epoch and fre- 
quency. The data between 1949 and 1976 in the frequency 
range 81.5 MHz to 9.4 GHz indicated a secular decrease of 
flux density estimated by Baars et al. (1977) to be frequency 
dependent and given by the relation 

4^ = 0.97(±0.04) - Q.m±Qm)log{v/GHz)% yr"\ (7) 
bat 

The most accurate data currently available covers the 
frequency range 38 MHz to 33 GHz and the period 1949 to 
2004 is given in Table [3] 

Fig. |6] shows the data from Table |3] divided into two 
epochs. The filled circles show the secular decrease from 
~ 1950 to 1970 while the crosses show the decrease from 
~ 1970 to 1990/2000. It is clearly evident that the earher 
data shows higher (~ 50%) secular decrease than the later 
data. The data from 1970 to 1990/2000 indicate a secular 
decrease of 



= °-68(±0.04) 



0.15(±0.04)log(!//GHz)% yr" 



(8) 



The material available is not sufficient to show whether 
the fall in the secular rate is linear with time or is or not 
(as proposed for the 81.5 MHz data by Hook et al. (1992)). 
Apart from the current (2001-2004) 33 GHz result the values 
given by equation 6 may be an upper limit to the secular 
variation of Cas A at the present epoch. The secular de- 
crease of Cas A will be discussed further in a related paper 
(Lancaster et al., in prep.). 

A frequency-dependent secular variation has implica- 
tions at the present epoch for the spectrum of Cas A. We 
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Figure 6. The secular decrease in the flux density of Cas A as 
a function of frequency for two epochs, 1950-1970 (filled circles) 
and 1965-2000 (crosses). The lines are least-square fits to the data 
at the two epochs. 



can apply equation (7) to the higher weight data from Table 
2 of Baars et al. (1977) and more recent data particularly at 
higher frequencies (Reichart & Stephens 2000, Mezger et al. 
1986, Mason et al. 1999, Liszt & Lucas 1999, Wright et al. 
1999, present paper) to derive spectra for 1965 (the Baars 
et al., epoch) and for 2000 at frequencies from 1 GHz to 250 
GHz). The decrease of the secular term with frequency as 
shown by equation (7) gives a flattening of the spectral in- 
dex with time. A synchrotron spectral index which flattens 
with time implies that particle acceleration is continuing in 
Cas A, presumably at the shock interfaces in the object (see 
e.g. Reynolds & Ellison 1992). 

4.3 Taurus A 

The Crab Nebula is the fiUed-centre remnant of the super- 
nova of 1054 AD. The central neutron star remains active 
and replenishes the supply of relativistic electrons in the 
nebula. There is some evidence for a small secular decrease 
of the fiux density at both radio and optical wavelengths. 
The present observations lead to a flux density of Tau A at 
33 GHz. Since the source is strongly linearly polarized we 
must determine the polarization in order to use Tau A as a 
calibrator at 1% accuracy. 

4-3.1 The spectrum of Tau A 

Baars et al. (1977) found a constant spectral index over the 
frequency range 1-35 GHz with value of a = —0.299 ±0.009. 
More recent data, including higher frequencies up to 350 
GHz are shown in Fig. [T] A constant spectral index of 
-0.278±0.007 flts the data between 1 and 350 GHz. The ex- 
pected turn-down of the spectrum is beyond this frequency 
range and is believed to be at 10 THz (Marsden et al. 1984). 
There is, however, some evidence for a small curvature of the 
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Figure 7. The log-log flux density spectrum of Tau A from 10 
MHz to 1000 GHz. A linear slope of a = -0.278 ± 0.007 fits the 
data from 1 to 350 GHz (dashed line). The slope between 1.8 and 
33 GHz (full line) is a = -0.351 ± 0.016. 



spectrum. We find a = -0.258 ±0.028 between 22 and 1800 
MHz, Q = -0.351 ± 0.016 between 1.8 and 33 GHz and 
a = -0.35 ± 0.03 between 8 and 347 GHz. Green (2002) 
gives the spectral index between 1.5 and 347 GHz of indi- 
vidual structures in the nebula to lie in the range —0.31 and 
—0.34. Similarly Mezger et al. (1986) give a spectral index 
between images at 10.7 and 250 GHz of a = -0.30 ± 0.05. 
We adopt a — —0.32 ±0.01 as an appropriate spectral index 
at « 30 GHz, which is consistent with the recent analysis 
by Macias- Perez et al. (2008). 



4.3.2 The VSA observations of Tau A 

Care is required when using Tau A as a calibrator in two 
regards. Firstly, it is an extended source with dimensions 
of 4 X 6 arcmin'^ at position angle (pa) ~ 140° and conse- 
quently requires a resolution correction when observed with 
the VSA in both the extended and compact configurations. 
Secondly, it is some 10% linearly polarized at frequencies 
of 30 GHz and above and this must be taken into account 
when observing with linear polarization. 

The observed optical polarization of Tau A is 9.2% at 
pa=159.6° (Oort and Walraven 1956). The integrated ra- 
dio polarisation increases from 2% at 1.4 GHz to 7% at 5 
GHz (Gardner & Davies 1966). The mean Faraday rotation 
measure (RM) in Tau A is —25 rad m~^ with values up to 
300 rad m~^ in some places. Any RM in this range leads to 
negligible Faraday rotation (FR) 30 GHz (a wavelength of 
1.0 cm) corresponding to a FR^ 1° (FR(X 1/A^). 

We now consider the VSA data for Tau A. 

(i) The linear polarization of Tau A. This can be esti- 
mated from long HA tracks of the source, remembering that 
the VSA is hnearly polarized at pa=0°. The maximum HA 
track for Tau A runs from HA=-3h 00m to +3h 00m. The 
observed data from the extended array on 26 December 2001 
calibrated against Gas A (and Jupiter) is shown in Fig.|8]by 
crosses. This data set is first corrected (circles) for resolu- 
tion effects of the VSA as for Gas A (section 4.2.2) using a 
5 GHz map of Tau; this model was confirmed by showing 
that the corrected visibilities were constant with baseline 
for both the compact and extended VSA configurations in 
the manner illustrated for Gas A in Fig. 5. The elevation 
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Figure 8. The flux density of Tau A plotted against HA for 
the observation on 26 December 2001 with the extended array. 
Crosses show the observed data calibrated against Jupiter; open 
circles show the data then corrected for resolution affects; flUed 
circles show the data further corrected for elevation (atmospheric 
absorption and the AGC effect). The resultant sloping HA plot 
is due to the linear polarization of Tau A. 



correction for absorption in the atmosphere and the AGC 
effect on the receivers (0.26 cosecB - see Section 3.4) were 
then applied giving the true variation of flux density with 
HA shown as filled circles in Fig. [8] The linearly polarized 
fiux density P is then estimated by fitting the corrected data 
to P cosec (2(p — 9)), where p is the polarization pa mea- 
sured east from north and 9 is the parallactic angle at the 
HA of observation. Data from 4 nights (26 December 2001, 
09 February 2002, 12 October 2002 and 25 January 2003) 
of good weather and analysed as above were averaged to 
determine the polarization of Tau A. The results of this fit- 
ting gave P = 24.6±1.6 Jy corresponding to 7.8±0.6% linear 
polarisation at pa=148°±3°. The errors quoted include the 
uncertainty in the atmospheric correction. 

(u) The fiux density of Tau A at 33 GHz. The total fiux 
density of Tau A is determined by the fitting procedure given 
in the previous section. The average of the 4 HA scans gives 
S(Tau A) = 322 ±4 Jy. The fiux density of Tau A at HA 
= Ohr (pa = 0°) is 332 ±4 Jy. A more extensive data set 
using short calibration tracks running from 3 January 2002 
to 11 September 2004 with the extended array gave a more 
accurate estimate. The individual data points were corrected 
for resolution and for the atmosphere and are plotted in 
Fig. |9] The use of Jupiter as the reference for Tau A led to 
the very low scatter (rms = 0.4%) in the data since they 
were both night-time objects and separated in time by less 
than 5 hrs. This data set gives a flux density for Tau A at 
pa = 0° of 332.8± 0.7 Jy on 1 January 2001. The secular 
decrease of Tau A indicated by Fig. is discussed in the 
next section. 



4.3.3 Secular decay of Tau A 

The 950-year old nebula is clearly in the decay phase al- 
though centrally reenergised from the spindown of its cen- 
tral pulsar. Theoretical models predict a secular decrease of 
0.16 to 0.4% per year at radio wavelengths (Reynolds and 
Chevalier 1984). An extended series of observations at 8.0 
GHz over the period 1969-1985 by Aller & Reynolds (1985) 
gave a decrease of 0.167±0.015% per year, a value consistent 
with the theoretical models (an interesting determination of 
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Figure 9. The flux density of Tau A measured relative to Jupiter 
for the period 3 January 2002 to 11 September 2004. The dashed 
line is a least-square fit to the data. The observed secular change 
is -0.22 ±0.07% yr-i. 



the optical decay of the continuum at 5000 Angstrom using 
data from 1960 to 2002.2 has been made by Smith (2003) 
who found a decay of 0.5±0.2% per year). 

The current 33 GHz observations of Tau A presented in 
Fig. 1^ and described in the previous section give a secular 
decrease over the period January 2002 to September 2004 of 

dS 



= -0.22 ± 0.07% yr'^. 
bat 



(9) 



This is very similar to the value obtained by AUer & 
Reynolds (1985) at 8.0 GHz for the longer but earlier pe- 
riod 1968 to 1985 and is also consistent with the models of 
Reynolds & Chevalier (1984). A secular decrease of ~ 0.2 % 
per year would seem to be appropriate for the present epoch 
at 33 GHz. 



4.4 NGC7027 

NGC7027 is young planetary nebula which is optically thick 
in Lyman continuum and is therefore an ionization-bounded 
ionized hydrogen region. The radio continuum and optical 
nebula have dimensions of 12 x 8 arcsec'^ (Basart & Daub 
1987). It is surrounded by a ring of neutral molecular gas 
and dust some 40 arcsec in diameter possibly extending to 
70 arcsec (see e.g. Hoare et al. 1992; Bieging et al. 1991). 
The expansion age of NGC7027 is 1200 yr with an expansion 
rate of 4.2±0.6 mas per year (Masson 1989). 

Three factors are of particular interest in the present 
study. The first is the activity within the nebula; its expan- 
sion velocity is 17.5±1.5 kms~^ with local gas velocities up 
to 55 kms"^ (Lopez et al. 1998, Cox et al. 2002). These 
motions may lead to secular variations of radio flux den- 
sity. The second effect is the dumpiness of the ionized gas 
in NGC7027 which results in optically thick knots up to 
at least 5 GHz (Bains et al. 2003). This determines the fre- 
quency range for which an optically thin radio thermal spec- 
trum (at Te ^14000 K) may be adopted. Thirdly, the dust 
within and surrounding NGC7027 may produce anomalous 
radio emission as in the case of the Helix Nebula (Casas- 
sus et al. 2004) and LDN 1622 (Casassus et al. 2006). Such 
emission, which peaks in the range 10-30 GHz (Draine & 
Lazarian 1998), needs to be considered when discussing the 
spectrum of NGC7027. 




10,0 100,0 
Frequency (GHz) 



Figure 10. The flux density spectrum of NGC7027 covering 
the range 408 MHz to 16000 GHz using data from the literature. 
The solid line is a fit consisting of a optically thick and thin free- 
free spectrum in addition to a single component modified black- 
body dust model. At frequencies of ~ 30 GHz, the spectrum is 
dominated by optically thin free-free emission. 



44.1 Spectrum of NG CI 021 

The radio spectrum of the integrated emission from 
NGC7027 is shown in Fig. [10] from the radio to the far in- 
frared (FIR). It consists of several components. At radio 
frequencies the spectrum is due to free-free emission, opti- 
cally thick at the lower frequencies and optically thin at the 
higher radio frequencies. The transition to optically thin oc- 
curs at ~5-10 GHz in agreement with the observation by 
Bryce et al. (1997) that the bright NW knot has an optical 
depth of ~2.0 at 5 GHz. A fit to the data between 10 and 80 
GHz gives an optically thin spectral index of — 0.119±0.025. 
The expected local value of the thermal spectral index for 
the observed electron temperature (14,000 K) in this fre- 
quency range lies between —0.111 (at 10 GHz) and —0.133 
(43 GHz) (Dickinson et al. 2003). 

The FIR emission consists of two components. One is 
warm (100 K) dust inside the 12x8 arcsec'^ ionized hydrogen 
region and the other is cool (30 K) dust in a neutral gas ring 
40 arcsec in diameter surrounding the nebula (Hoare et al. 
1992, Bieging et al. 1991). Both these dust components are 
covered by the observing beamwidth of the VSA and are 
included in the integrated emission described here. 

The question of whether there is significant anoma- 
lous radio emission from either or both of the dust compo- 
nents can be addressed. There is no obvious emission with 
a peaked spectrum of roughly an octave width centred at 
~ 20 GHz in the spectrum of Fig. [TD] For a normal cool dust 
cloud the expected emission based on the lOO^tm (3 THz) 
flux density would be « 0.3 Jy (Davies et al. 2006, Wat- 
son et al. 2005, Casassus et al. 2006). The emission from a 
normal HII region (Dickinson et al. 2006 & 2007) would be 
« 0.1 Jy. The latter is more consistent with Fig. [TOl 
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Figure 11. The flux density of NGC7027 calibrated by Jupiter 
and Cyg A for the period June 2002 to July 2003, as measured by 
the VSA. There is weak evidence for a higher flux density in the 
first 6 data points. The annual change for the full data set (dashed 
line) is —0.8 ± 0.3% yr~^\ the change for the data excluding the 
first 6 points (full Unc) is -0.31 ± 0.37% yr"!. 
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Figure 12. The flux density of NGC7027 from published data 
for two frequency ranges over the period 1967 to 2004 with least- 
square fits. 

(a) Full line and filled circles; optically thin frequencies, 8.0 to 
43 GHz, and 

(b) dashed line and crosses; partially optically thick frequencies, 
4.8 to 6.6 GHz. 



4.4.2 The variability of NGC7027 

The flux densities measured over the period June 2002 to 
July 2003 are shown in Fig. 1111 There is weak evidence for 
a higher flux density in the 6 observations before 10 July 
2002. After this time the flux density remains constant for 
the year; a linear flt to the latter data gives an insignificant 
annual change of —0.31 ±0.37% yr~^. A flt to the full run of 
the data gives a decrease of 0.8±0.3% per year. The average 
flux density of NGC7027 in the period June 2002 to July 
2003 is 5.39±0.04 Jy. 

The long-term secular variation of NGC7027 at ra- 
dio frequencies is often discussed (see e.g. Ott et al. 1994; 
Masson 1989). A substantial data set covering the fre- 
quency range 4.8 to 43 GHz is available in the literature 
extending back to 1967. Fig. [T^ shows the flux density 
of NGC7027 from 1967 to 2003 in two frequency ranges: 
4.8 — 6.6 GHz where free-free self-absorption becomes sig- 
nificant and 8.0 — 43 GHz where the spectrum is optically 
thin as can be seen in Fig. 1101 All flux densities are cor- 
rected to a frequency of 33 GHz, assuming an optically thin 
spectral index a = —0.12 as expected theoretically for this 
frequency range (Dickinson et al. 2003). The 4.8 - 6.6 GHz 
flux densities (median frequency 5.2 GHz) are 82 ± 3% of 
the optically thin frequencies, indicating r=0.18 at 5 GHz. 
Giving each point equal weight, the secular change in the 
flux density at optically thin frequencies (8.0 to 43 GHz) is 
J c 

= -0.159 ± 0.043% yr-\ (10) 

bat 

At the slightly optically thick frequencies (4.8 - 6.6 
GHz) 

= -0.45 ±0.13% 2/r--i (11) 

bat 

when taking all the data. However, discarding the ear- 
lier data with larger errors, the secular change from 1975 to 
1996 is 

4^ = -0.183 ± 0.052% yr'^ (12) 
bat 

in close agreement with the optically thin data. We 
therefore have consistent results from two independent data 



sets. On combining the above data, we conclude that for 
NGC7027 

4^ = -0.17 ± 0.03% (13) 
bat 

over the period ^^1970 to the present. 

After completion of this work Zijlstra, van Hoof & Per- 
ley (2008) published observations made with the VLA at fre- 
quencies in the range 1.275 to 43.34 GHz covering the period 
1983-2006. They find that over this shorter period the flux 
density increases at lower frequencies where the emission is 
optically thick and decreases at higher (optically thin) fre- 
quencies. The turnover occurs between 2 and 4 GHz. Their 
value of decrease at optically thin frequencies of 0.145 ± 
0.005 % yr~^ is consistent to that of the present study cov- 
ering data from the 36 year period to 2003. 

4.5 Venus 

The radio spectrum of Venus is the result of contributions 
from different depths in its atmosphere at high frequencies 
and from the surface at low frequencies. At the surface the 
temperature is « 750 K and the pressure ~ 90 bar. The 
millimetric emission arises from the higher atmosphere at 
h~ 80 km where the effective temperature is 300 K. At inter- 
mediate radio frequencies the emission arises from the finite 
microwave opacity of SO2 and H2SO4. At 33 GHz where the 
brightness temperature is « 460 K the effective altitude of 
the emission is 35 km; this is the altitude at which space 
probes recorded a temperature of 460 K and a pressure of 
« 5 bar. At altitudes of ~100 km a strong Venusian di- 
urnal effect has been measured by space probes. We have 
examined our data to determine whether this diurnal effect 
extends down to altitudes around 35 km where the 33 GHz 
emission arises. 



4.5.1 The spectrum of Venus 

Fig. I13f a) shows the radio brightness temperature spectrum 
of Venus. The low frequencies come from the surface region 
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Figure 13. (a) The brightness temperature spectrum of Venus 
from the literature. The full line shows the least-square best fit 
slope between 10 and 100 GHz which gives a T;, spectral index 
at 33 GHz of -0.278 ± 0.026. 

(b) The brightness temperature spectrum of Saturn from the lit- 
erature. 



where the temperature is ^^750 K. The higher GHz frequen- 
cies arise in successively higher levels in the Venusian atmo- 
sphere where temperatures are lower. The slope of the T;, 
spectrum at 33 GHz, estimated from a best fit to the range 
10 to 100 GHz, is -0.278 ± 0.026 (the temperature spectral 
index, /3, is related to the fiux density spectral index, a, by 
P = a-2). 



4-5.2 Brightness temperature variation with illumination 

Fig. [14] shows the 33 GHz brightness temperature of Venus 
as a function of illumination for the period August 2002 to 
September 2004. During this period, Venus went through 1.5 
Synodic cycles. Data are shown separately for eastern and 
western illumination, corresponding to times when the Venu- 
sian surface was entering or leaving the long period of dark- 
ness. No difference is seen between these two phases of illu- 
mination of more than 2 K ( 0.3 %). The variation of bright- 
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Figure 14. The brightness temperature of Venus plotted as func- 
tion of illumination for the period August 2002 to September 
2004. Filled circles show eastern illumination (surface entering 
darkness) and open circles show western illumination (leaving 
darkness). The full line is the best linear fit to the data corre- 
sponding to +1.5 ± 1.4% change from zero to full illumination. 



ness temperature with illumination suggests a weak increase 
of temperature with illumination of -f 1.5±1.3% (6.9±6.0 K) 
from zero to full illumination. 



4.6 Saturn 

Saturn is used as a secondary calibrator by the VSA. For 
accurate measurements the effect of the rings must be taken 
into account. The rings themselves emit with a brightness 
temperature of 10-20 K and also block out the emission from 
the disk of Saturn (see for example Janssen & Olsen 1978, 
Conway 1980; Ulich 1981). As a consequence, the effective 
brightness temperature of Saturn varies with the tilt angle 
B (the Saturnicentric latitude of the earth). We present data 
for the period June 2002 to August 2003. 



4-6.1 The spectrum of Saturn 

The microwave spectrum of Saturn taken from the literature 
is plotted in Fig. I13r b) in terms of its brightness tempera- 
ture. There is some uncertainty as to the contribution of 
the rings to the data presented. Since the ring effect is typi- 
cally less than 5 percent, it lies within the error assigned to 
most data points in Fig I13f b) . The shape of the spectrum 
is quantitatively similar to that of Jupiter - increasing from 
~ 130 K at mm wavelengths to 540 K at 94 cm. As in the 
case of Jupiter there is evidence for the presence of ammonia 
as the principle source of opacity at radio frequencies with 
its inversion band at « 24 GHz. At 30 GHz the frequency of 
interest in the present study, the effective temperature spec- 
tral index is /? = 0.00±0.05. At frequencies of 10 - 43 GHz, 
/3 = -0.07±0.05, while at 10-200 GHz, /3 = -f0.041±0.015. 



4-6.2 Secular variation of Saturn at 33 GHz 

The VSA programme gave an accurate data set for the inte- 
grated flux density of Saturn extending from 29 June 2002 
to 02 August 2003. During this period the tilt angle moved 
between B = 26° and 27° with a quasi - annual period. 
The brightness temperature is then calculated by assuming 
the observed flux density arises from Saturn's disk whose 
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Figure 15. The brightness temperature of Saturn at 33 GHz for 
the period 29 June 2002 to 02 August 2003. There is no secular 
variation greater than 0.4% yr~^ over this period. 



dimensions are given in the Ephemeris. Fig [15] is the bright- 
ness temperature plot over this period and shows no secular 
variation greater than 0.4% per year at 95% confidence level. 
Our results indicate that the mean brightness temperature 
of Saturn (disk plus rings) between June 2002 and August 
2003 is 140.50±0.12 K; for the data set the mean tilt angle 
B = 26.5°. 



5 CONCLUSION 

We have determined the flux density or brightness tempera- 
ture at 33 GHz of 7 sources used by the VSA for calibration; 
the results are summarised in Tables|4]&[5l This list includes 
sources widely used for calibration in this frequency range in 
radio-source and CMB studies; they are the strong sources 
and planets. The high level of reliability of the VSA and the 
stable atmosphere at the Izana site have assured accurate 
relative measurements of these sources. We use Jupiter as 
the reference standard in this work and adopt a brightness 
temperature of 146.6±0.75 K at 33 GHz as determined in 
the WMAP survey (HiU et al. 2008). Any refinement in the 
future can then be applied to the present results. In order to 
make our results useful for accurate flux density/brightness 
temperature measurements in the frequency band around 33 
GHz we derive a spectral index from published data. 

Since we have observations of the calibrators extending 
over nearly 4 years, we are able to make an accurate assess- 
ment of the calibrator variation over this period. With rms 
individual comparison measurements of 0.4 to 1.0 % depend- 
ing on the time separation of the pair of observations, we are 
sensitive to small secular flux density or brightness tempera- 
ture variations. A significant secular decrease in fiux density 
was established for Gas A and Tau A over the 4-year observ- 
ing period. By combining our data with historical data in 
the literature, we have also demonstrated a significant secu- 
lar decrease for NGG7027 at optically thin frequencies. Over 
this period the secular variation of the two sources expected 
not to vary, Gyg A and Jupiter, was less than 0.1% yr"'^. 
Venus showed a slight (1.5±1.3%) variation with illumina- 
tion. 
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Radio source calibration at 30 GHz 13 



Table 4. Secular change in flux density (or brightness temperature) of calibrators relative to Jupiter at 33 GHz (epoch 2001.0). 



Radio sources 


dS/Sdt (% yr-i) 


Gas A 


-0.394 ±0.019 


Tau A 


-0.22 ±0.07 


NGG70271 


-0.16 ±0.04 


Cyg A 


+0.043 ± 0.039 


Venus ^ 


+1.5 ± 1.3 



1 At optically thin frequencies; see Section 4.4.2. 

2 Secular variation of Tj, with illumination; see Section 4.5.2. 



Table 5. Summary of calibrator flux densities or brightness temperatures at 33 GHz and epoch 2001.0, all relative to an adopted Jupiter 
= 146.6 K. 



Radio sources Flux density Error Spectral index 

(Jy) (Jy) (at 33 GHz) 

Gas A 182.0 0.1 -0.69 + 0.02 

Gyg A 36.4 0.2 -1.208 ±0.017 

Tau Ai'2 322 4 -0.32 ± 0.01 

NGG70273 5.39 0.04 -0.119 ±0.025 



Planets Brightness temperature Error Spectral index 

(K) (K) (at 33 GHz) 

Jupiter 146.6 0.75 +0.248 + 0.004 

Venus 460.3 3.2 -0.278 + 0.026 

Saturn"' 140.50 0.12 0.00 + 0.05 

^Tau A is 7.8+0.6% yr'^ Unearly polarized at pa = 148° + 3°. 
^Tau A flux density is 332.8+0.7 Jy at pa = 0°. 
^Flux density at 2003.0. 

"Flux density at 2003.0; tilt angle B of rings is 26.5°. 
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